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Abstract 

The ^‘^Mg(p, Q;)^^Na reaction was measured at the Holifield Radioactive Ion Beam Facility at Oak 
Ridge National Laboratory in order to better constrain spins and parities of energy levels in ^^Na 
for the astrophysically important ^'^F(a,p)^‘^Ne reaction rate calculation. 31 MeV proton beams 
from the 25-MV tandem accelerator and enriched ^^Mg solid targets were used. Recoiling ^He 
particles from the ^^Mg(p, Q;)^^Na reaction were detected by a highly segmented silicon detector 
array which measured the yields of ^He particles over a range of angles simultaneously. A new 
level at 6661 ± 5 keV was observed in the present work. The extracted angular distributions for 
the first four levels of ^^Na and Distorted Wave Born Approximation (DWBA) calculations were 
compared to verify and extract angular momentum transfer. 
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I. INTRODUCTION 


Detecting 7 -rays from the decay of the long-lived radionuclide "^^Ti ^i/ 2 = 59 . 1 y) provides 
a direct calibration of the nucleosynthesis in core-collapse supernovae jj . Because of its im¬ 
portance, the abundance of ^“^Ti has been the focus of many studies jJ?]. As an example, 


G. Magkotsios et al. investigated the "^"^Ti abundance produced from core-collapse super¬ 
novae by studying the impact on the ^^Ti abundance evolution of variation of the nuclear 
reactions, including ( 0 , 7 ), (a,p), (p, 7 ), (p, a), {Py^), and (a,n) in light and intermediate - 
mass targets. In their sensitivity study, it was found that the variation in the ^^F(a,p)^°Ne 
reaction rate causes a “primary” impact on the ^^Ti abundance. The ^^F(a,p)^°Ne reaction 
rate, however, has never been measured. Because the reaction rate may be dominated by 
the properties of energy levels of ^^Na above the a-threshold at 6.561 MeV, searching for 
energy levels of ^^Na and studying their properties may impact our understanding of the 
abundance evolution of ^^Ti. 

The ^^Mg(p, a)^^Na reaction was measured at the Holiheld Radioactive Ion Beam Facility 
(HRIBF) at Oak Ridge National Laboratory (ORNL) in order to make a spectroscopic study 
of the energy levels in the ^^Na for the ^’^F(Q;,p)^°Ne reaction rate at stellar temperatures. 
The ^^Mg(p, a)^^Na reaction was reported only once previously, by J. G. Pronko et al. j^, 
where the lower excited states < 5 MeV) of ^^Na were studied. Spins and parities of 
three excited states located at 0.332, 1.716, and 2.829 MeV were constrained by analyzing 
particle- 7 -ray angular correlations. In the present work, energy levels of ^^Na up to ~ 8.5 
MeV were observed, which will provide the hrst spectroscopic results of the ^"^Mg(p, Q;)^^Na 
reaction in the energy range of E^ = 5-8.5 MeV. An partial analysis of the data, for the 
lowest four states, is presented in this contribution. 


II. EXPERIMENT 

A schematic diagram of the experimental setup is shown in Fig. [H A 31 MeV proton 
beam from the 25-MV tandem accelerator bombarded a pure ^“^Mg solid target placed at 
one of two different target positions at 3.75 in. and 1.75 in. from the front side of the 
silicon detector array to subtend two separate but overlapping angular ranges. The targets 
were isotopically enriched (99.9 %) and had thicknesses of 520-/Ug/cm^ and 516-/Ug/cm^, 
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FIG. 1: A schematic diagram of the experimental setup is shown. A 31 MeV proton beam impinged 
on an enriched ^“^Mg solid target. Recoiling a particles from the ^^Mg(p, a)^^Na reaction were 
detected at forward angles by the SIDAR. A thick graphite beam stop was used for beam current 
monitoring. 

respectively. A 0.375 in. thick aluminum plate with an aperture (~ 0.75 in. diameter) was 
placed just upstream of the target ladder in order to prevent the fragile silicon detectors 
from being exposed to the intense proton beam. 

Alpha particles from the ^‘^Mg(p, a)^^Na reaction were detected by a large area annular 
silicon detector array (SIDAR) 9[. The SIDAR was composed of four trapezoidally shaped 
AE — E telescopes for particle identihcation. Each telescope was conhgured with a thin (100- 
pm) AE detector backed by a thick (1000-/im) E detector. With this detector thickness 
and a beam energy of 31 MeV, even the most energetic a particles from the (p, a) reactions 
(i.e., the a particles associated with the ground state of ^^Na) can be fully stopped in the 
detector system. The detector wedges were tilted upstream 43° from the perpendicular to the 
beam axis in order to cover a larger angular range. Since each detector is segmented into 16 
azimuthal strips, the reaction cross sections were measured at multiple angles simultaneously, 
and thus the angular distributions of recoiling a particles were extracted without the need 
for relative normalizations. The angles covered by the SIDAR were 17° < 9iab < 44° (19° < 
9c.m. < 50°) and 24° < 9iab < 63° (27° < 9c,m. < 70°) for the two target positions. By placing 
two targets at different positions, a wider range of angles could be covered without changing 
any other parameters of the experimental setup. Three out of four detector telescopes were 
with respect to gain for the ^^Mg(p, (i)^^Mg reaction which was measured simultaneously. 


3 






























100 200 300 400 500 

AE + E (channel) 

FIG. 2: (Color online) An example of the energy loss (AFl) versus total energy {/S.E + E) plot 
obtained at 6iab = 27.2° {9c.m. = 30.4°) is shown. Events falling in the red gate were identified as 
a particles. Other reaction products are also labeled. 

Since the energies of the a particles from ^‘^Mg(p, Q;)^^Na are abont factor of 2 higher than 
those of the denterons from the ^^Mg(p, d)^^Mg reaction, the fourth telescope was optimized 
in gain for the (p, a) channel. The beam current was continuously integrated from a thick 
graphite beam stop which was placed on the downstream side of the target chamber. 

III. DATA ANALYSIS 


To calibrate the energy response of the silicon detectors, an a emitting source composed 
of ^^'^Pu (5.157 MeV), ^^^Am (5.486 MeV), and (5.805 MeV) was used. By using a 

peaks of the three energies, the energy offset of each Analog to Digital Converter (ADC) 
channel could be determined as well. An additional calibrated ^^"^Cm source was used to 
measure the solid angle subtended by each strip. The measured solid angles were cross¬ 
checked with calculated solid angles using the known detector geometry. The measured 
solid angles agreed with geometric calculations to within 3 %. 

Light charged particles were identihed by a standard energy loss technique. A typical 
particle-identihcation plot is shown in Fig. [2l The plot was obtained for the strip subtending 
at 6iab = 27.2° {9c.m. = 30.4°) in this case. Events falling in the red gate in Fig. 2 were 
identihed as ^He particles. The ^He yields were clearly identihed as shown in the hgure 
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FIG. 3: The number of counts per channel versus a energy plot at Qiah = 27.2° is shown. Excitation 
energies are labeled in MeV. The labeled levels have well-known spin assignments. 


without significant evidence of contamination from the other charged particle groups such 
as p, d, t, and ^He. The total energies of the ^He particles were reconstructed by summing 
the energies deposited in the AE and E detectors. Fig. [3] shows an example of a total 
energy spectrum of gated “^He particles from the particle identihcation plot. Many energy 
levels of the radioactive ^^Na nucleus are evident in this spectrum. 

Because the recorded ADC channels are not always linearly dependent on the actual 
energies of incident particles, internal energy calibrations are necessary to more precisely 
determine the energy spectra. The internal calibrations were performed at each angle using 
two energy levels of ^^Na located at 0.332 MeV and 2.829 MeV excitation, both of which 
were strongly populated as shown in Fig. [3l Calibrated ^He energies were then converted to 
excitation energies in ^^Na using the well-known detector geometry and reaction kinematics. 
Five levels, including the ground state, were clearly identihed as labeled in Fig. |3l As shown 
in the figure, however, many excited states were populated as well at the energies above 
Ex > 2.829 MeV. Improved internal energy calibrations are required to determine precise 
excitation energies of the higher-lying states due to the required extrapolation to low ^He 
energies. 

In order to verify the results of the internal calibrations, a newly developed relativistic 


0 . 


By using the 


kinematics calculation code VisKin (Visualized Kinematics) was used 
object-oriented computer programming language Java |^, VisKin helps users to interact 
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FIG. 4; (Color online) The result of internal calibrations performed at Oi^b = 


29.1° data are shown. 


with the calculated results so that the internal calibrations can be performed relatively 
easily. Upon execution, VisKin requires users to upload an image of an energy spectrum 
that will be used for the internal calibration. Information regarding the nuclear reaction such 
as masses and atomic numbers of particles participating in the reaction, beam energy, and 
emitting laboratory angle of the ejectile is also required to perform the relativistic kinematics 
calculations. The program then reads in the excitation energies of the heavy recoil, which 
are extracted from the Evaluated Nuclear Structure Data File (ENSDF) database provided 
by the National Nuclear Data Center of the Brookhaven National Laboratory 13|. After 
calculating the energies of the ejectiles that would be produced by the population of known 
energy levels of the heavy recoils, the program displays the energies as red lines on top of 
the energy spectrum image which the user provided as an input. The calculated results (red 
lines) can be translated and zoomed so that the internal calibrations can be determined 
easily (Fig. H]). 

The result of internal calibrations performed for the ^^Mg(p, Q;)^^Na reaction study is 
shown in Fig. 01 The experimental data represent counts per channel obtained at 9iab = 
29.1° as a function of the total a energy. As shown in the hgure, lower-lying levels of ^^Na 
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FIG. 5: (Color online) Excitation energies extracted at 11 angles for a peak located near channel 
number 2770 in Fig. 0]are plotted as a function of Oc.m.- The black squares (red triangles) represent 
the calculated excitation energies of ^^Na (^^N) assuming the peak arose from the ^‘^Mg(p, a!)^^Na 
(^®0(p, Q!)^^N) reaction. 


up to Ex ~ 6.9 MeV were well reproduced by known energy levels from references [8|,|]J4l7| 
except the level located near channel number 2770. The excitation energy in ^^Na for this 
level was calculated to be 6661 ± 5 keV, which does not correspond to any known 


levels. The excitation energy uncertainty was calculated by using Eq. (1) of Ref. jl8 |. 
Excitation energies extracted at 11 angles are plotted in Fig. |5] as a function of 9c.m- As 
shown in the figure, the excitation energies are rather consistent indicating that the peak 
did not arise from any contamination in the target. 

The differential cross sections of the identified levels in the center of mass system at each 
angle were calculated as 


(E] ^ ( 1 ) 

\da)., INMh' '■ ’ 

where Yj^o is the yield of a particles from the ^^Mg(p, a)^^Na reaction for an energy level 
2 emitted at an angle of Q in the center of mass system, I is the number of beam particles 
incident on the target, N is the number of target atoms per unit area, and AVlg is the solid 
angle covered by a strip of SIDAR in the center of mass system. 

Angular distributions of four levels located at energies of 0.0, 0.332, 1.716, and 2.424 
MeV were extracted from the experimental data. Since these four levels have well-known 
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TABLE I: The Optical model parameters used for DWUCK4 code in this work are shown. The 
definitions of these parameters follow the conventions in the Ref. 0. 


Particle VR (MeV) 

^OR (fm) 

AR (fm) 

VI (MeV) 

ro/ (fm) 

AI (fm) 

roc (fm) 

P 

44.53 

1.141 

0.15 

17.51 

1.26 

0.64 

1.3 

^He 

90.1 

1.2 

0.87 

13.8 

1.8 

0.99 

1.3 



FIG. 6: (Color online) The angular distributions of a particles from the ^^Mg(p, Q:)^^Na reaction 
(circles) and best-htting DWBA calculations (black solid lines) for four levels are shown. The 
statistical uncertainties of differential cross sections are smaller than the size of the circles. All 
excitation energies, labeled on the bottom-left of each plot, are in MeV. 


spin and parity assignments 


14 


ITj , comparing experimental angular distributions of the 


levels with Distorted Wave Born Approximation (DWBA) calculations considering proper 
/ transfer will verify the validity of the optical model parameters used in the calculations. 
DWBA calculations were performed using the zero range computer code DWUCK4 19|]. The 
optical parameters were adopted from previous work 2^ and slightly modihed to better £t 
experimental data. If extracting spectroscopic factors is the main purpose of the study, 
hne-tuning optical model parameters obtained from previous scattering experiments would 
be inappropriate. In the present work, however, the goal is to find the proper I transfer that 
is consistent the maxima and minima observed in the experimental angular distribution. 
Table, m summarizes the optical model parameters used for the calculations. 

Fig. [6] shows the angular distributions of '^He particles for the first four levels including 
the ground state. The circles represent experimental data obtained from this work, and black 
solid lines represent DWBA calculations that best fit the experimental angular distributions. 




















Statistical uncertainties are marked as error bars. In most cases, however, the uncertainties 
are smaller than the size of data points. Triton transfer was assumed for the calculations. 
Transferred I values are also shown in each plot. The number of beam particles incident 
on the target (/) and the number of target atoms per unit area {N) in Eq. ([T]) are not 
included in the calculations of the experimentally obtained differential cross sections. This 
does not affect the shape of the angular distributions, since the cross sections were measured 
at all angles simultaneously. As shown in Fig. |6l the DWBA calculations well reproduce 
the extracted angular distributions for the four well-known energy levels. 


IV. CONCLUSION AND FUTURE PLAN 

A spectroscopic study of the ^'^Mg(p, a)^^Na triton transfer reaction was made by using 
31 MeV proton beams at the HRIBF of the ORNL to constrain spins of energy levels in 
^^Na for the astrophysically-important ^^F(q;, p)^°Ne reaction rate. The energies and angu¬ 
lar distributions of recoiling “^He particles from the reaction were measured using a highly 
segmented silicon detector array. By comparing experimental differential cross sections of 
strongly populated ^^Na energy levels with theoretical DWBA calculations, we could verify 
the following spin and parity assignments: ground state - | , i?3.=0.332 MeV - | , F'j.=1.716 
MeV - 1^, F'a,=2.424 MeV - which supports the validity of the optical model parameters 
used for the calculations. 

Many energy levels in ^^Na up to ~ 6.9 MeV were observed, some of which do not have 
dehnite spin assignments and others have not been observed previously. To determine precise 
excitation energies and to constrain spins and parities of the higher-lying states, improved 
internal calibrations are required. Analysis of the experimental data obtained at the second 
target position will extend the angular coverage up to Oc.m. ~ 70°, which will be useful to 
constrain the spin assignments. The ^^F(a:, p^^Ne reaction rate will be updated at stellar 
temperatures as well. 
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